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ABSTRACT 

We present a timing and broad-band pulse phase resolved spectral analysis 
of the transient Be X-ray binary pulsar 1A 1118-61 observed during its out- 
burst in January 2009 using Suzaku observations. Suzaku observations were 
made twice, once at the peak of the outburst, and once 13 days later at its 
declining phase. Pulse profiles from both observations exhibit strong energy 
dependence with several peaks at low energies and a single peak above ~10 
keV. A weak, narrow peak is detected at the main dip of the pulse profiles from 
both observations in the energy bands below 3 keV, indicating the presence of 
a phase dependent soft excess in the source continuum. The broad-band en- 
ergy spectrum of the pulsar could be fitted well with a partial covering cutoff 
power-law model and a narrow iron fluorescence line. We also detect a broad 
cyclotron feature at ^50 keV from both observations which is a feature com- 
mon for accretion powered pulsars with high magnetic field strengths. Pulse 
phase-resolved spectral analysis shows an increase in the absorption column 
density of the partial covering component, as well as variation in the covering 
fraction at the dips of the pulse profiles, that naturally explains the energy 
dependence of the same. The cyclotron line parameters also show significant 
variation with pulse phase with a ~ 10 keV variation in the cyclotron line 
energy and a variation in depth by a factor of three. This can be explained 
either as the effect of different viewing angles of the dipole field at different 
pulse phases, or a more complex underlying magnetic field geometry. 

Key words: X-rays: binaries- X-rays: individual: 1A 1118-61- stars: pulsars: 
general 
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1 INTRODUCTION 

Be/ X-ray binaries are a subtype of neutron star high mass X-ray binaries (HMXBs) where 
the optical counterpart is a dwarf, subgiant or a giant O/Be star (luminosity class III-V) 
Reigl (120111 ). These systems typically have a large orbital period, moderate eccentricity, 
and being a non-supergiant system, the Be star lies deep within the Roche lobe. X-ray 
emission is believed to be due to the accretion from the circumstellar disc of the rapidly 
rotating Be star onto the compact star. Most Be/ X-ray binaries are transient systems and 
display two types of outbursts such as (i) short, periodic, less luminous Type I outbursts 
(L x 10 36 — 10 37 erg s -1 ) lasting for a few days occurring at the phase of the periastron 
passage, and (ii) giant, longer, Type II outbursts (L x ^ 10 37 erg s _1 ) lasting for several 
weeks to months occurring when a large fraction of matter is accreted from the Be star's 
disk. Type II outbursts are sometimes followed by smaller recurrent type I outbursts. 

Almost all Be/ X-ray binaries are accretion powered pulsars having a wide range of 
periods from seconds to hundreds of seconds. These are highly magnetized neutron stars 
accreting matter, channeling it along magnetic fields onto the magnetic pol es, and forming 



an accretion column of X-ray em i tting plasma on the neutron star surface ( Pringle &: Rees 



1972 



Davidson fc Ostrikcr 



1973 



Gnedin fc Sunyaev 



1973 



Lamb. Pethick fc Pines 



19731). 



The spectrum of these o bjects can be described by Comptonization of the so 
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Brainerd fc Meszaros 
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t X-rays by scat- 



Meszaros fc Nagel 



b|). In the presence of such strong magnetic fields, photons undergo 
resonant scatterings by electrons in Landau levels to produce the Cyclotron Resonance Scat- 

These features have been found in more than 16 accretion powered 



tering F eatures (CRSF 
I 

1995 



pulsars ( IMihara 



Coburn. Heindl fc Rothschild 



20021 ) . The cyclotron line properties 



provide valuable information on the details of the emission, electron temperature, optical 
depth, and viewing geometries, and can map the magnetic field structure. For a brief re- 
view on the t empo ral and spectral properties of transient Be/X-ray binary pulsars, refer 
Paul fc Naik! f ]201lh . 

1A 1118-615 is a hard X-ray tra nsient pulsar that was discovered with Ariel V in an 
outburst in 1974 ( lEyles et al.lll975l ). The same series of observatio ns revealed X-ray pul- 



sations with a duration of 405.6 s ( Ives. Sanford fc Burnell 



19751 ) . The optical counter- 
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part (He 3-640= "WRA 793 ") is a highly reddened Be star classified as an 09.5IV-Ve star 
( iChevalier fc Ilovaiskylll975l). It shows strong Balmer emission lines indicating the presence 
of an extended envelope ( Motch et al. Il988 ). The distance to the source is estimated t o 
be 5 ± 2 kpc from photometric and sp ectroscopic observations ( Uanot-Pacheco et al.lll98ll ). 
UV observations ( ICoe fc Paynd Il985l ) reveal a P-Cygni profile in the C IV line which in- 
dicate s a stellar outflow with a velocity of the order of 1600±300 km s _1 . ICoe fc Payne 
(119851 ) also confirmed the magnitude of the companion and from the extinction value, the 
distance was suggeste d to be about 4 kpc. An orbital period of 24±0.4 days was reported by 



Staubert et al. 



(120111 ). There have been thr ee giant outbursts so far detected for this source. 
The first was in 1974 (IMarashi et al.l 119761 ) . The second was detected in January 1992, and 
was observed by the B urst and Transient Source Experimen t (BATSE) covering an energy 



range of 15 — 50 keV (ICoe at al. 



1994 



Marashi et al. 



19761 ). The outbur st lasted for ~3 



days and pulsations with a period ~406.5 s were detected up to 100 keV (ICoe at al.lll994j ). 
The same work also reported a spin change rate of -0.016 s/day at the decay of the outburst. 

The source remain ed in quiescence for 17 years until 4 January 2009 when a third outburst 
was detected by Swift (jManganoll2009l ). Pulsations were detected at ~ 407.68 s with the pulse 
profile showing a complex structure. The complete outburst was regularly monitored with 
RXTE (Rossi X-ray Timing Explor er). INTEGRAL de tected a flaring activity in the source 
~ 30 days after the main outburst (ILeyder et al.l 120091 ) . Suzaku observed the source twice, 
once at the peak of the outburs t and o nce 13 days later. RXTE observations of the source was 

(120101 ) who presented the broad-band spectrum, and reported 



analyzed by 



Doroshenko et al. 



the presence of a broad prominent absorption feature at ~55 keV which they interpreted as a 
CRSF. Detailed timin g and spectral analysis of the same set of observations was reported by 



Devasia et al. 



(l2011al ) which included the energy dependence of the pulse profiles and pulse 
phase-resolved spectroscopy of the continuum spectra. They also detected quasi-periodic 
oscillations (QPO) at 0.07-0.09 Hz which showed a significant energy dependence. Analysis 

(120 111 ) which confirmed the 



Suchy et al 



of the same Suzaku observations was reported by 
presence of the CRSF at ~55 keV. They also performed pulse phase-resolved spectroscopy 
for three different pulse phases i.e. the two peaks and the minimum of the pulse profile. 
Suchy et al.l (120111 ) thus reported a change in the continuum parameters like the photon 
index and high energy cutoff, and also the centroid energy in the three above mentioned 
phases. 

The pulse profiles of 1A 1118-615 show multiple energy dependent dips. Dips in the 
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pulse profiles can be naturally produced by absorption in the accretion stream that is 



phase locked with the n eutron star. This had been explore d for GX 1+4 (jGallowav et al. 



20011), RX J0812 4-3114 flCorbet fc Peelell2000l ). 1A 1 118-61 



1 (IDevasia et al. 



Devasia et al. 



201 lbh with RXTE, KS 1947 + 300 flNaik et al 



Suchy et al 



200i]) with Beppo-Sax, and 



2011ah . GX 304- 



( 1201 ll ) did include a phase de- 



GRO J1008-57 with Suzaku flNaik et al.ll201ll ). 
pendent additional absorption component in their analysis but did not probe the narrow dips 
of the pulse profiles, as they carried out phase reso lved analysis by dividing the pulse profile 
into only three phase bins. IDevasia et al.l ( l2011al ) used RXTE-PCA (RXTE Proportional 
Counter Array) data covering a limited energy-band. We report here a more comprehensive 
timing and broad-band spectral analysis of Suzaku observations of 1A 1118-61 during the 
2009 outburst. We have probed the energy dependence of the pulse profiles in more detail 
and have carried out pulse phase-resolved spectroscopy in narrow phase bins to explain its 
complex nature. 

We have also carried out pulse phase-resolved spectroscopy of the cyclotron absorption fea- 
ture in this source and obtained one of the most detailed pulse phase-resolved information 
available for CRSF till date. In addition, we have detected a narrow phase dependent soft 
excess peak in both observations of this source. In section 2, we describe details of the 
observations and data reduction. In section 3, we describe the timing analysis including 
the energy dependence of the pulse profiles and the power density spectra. In section 4, 
we present the pulse phase averaged and pulse phase-resolved spectroscopy of the Suzaku 
observations followed by discussions & conclusions in section 5. 



2 OBSERVATIONS & DATA REDUCTION 



The hard X-ray transient pulsar 1A 1118-61 was observed with Suzaku, once at the peak of 
the outburst and once in its decline 13 days later. Figure [T] shows the one day averaged long- 
term light curve of 1A 1118-61 using data from the the All Sky Monitor (ASM) (1.5 — 12 keV) 
onboard RXTE and from the Swift-B AT all sky monitor (15 — 50 keV) Q. The arrow marks 
i n the figure shows th e Suzaku observations of the pulsar during the 2009 outburst Suzaku 
( iMitsuda et al.l 120071 ). covers the 0.2-6 00 keV energy range . It has two sets of instruments, 
the X-ray Imaging Spectrometer XIS ( IKoyama et al.l 120071 ) covering the 0.2-12 keV range, 



1 See ftp:/ /heasarc. gsfc.nasa.gov/xte/data/archive/ASMProducts/ and http:/ /heasarc. nasa.gov/docs/swift/results/transients/ 
for details 
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the Hard X-ray Detector (HXD) having PIN diodes (ITakahashi et al.l 120071 ) covering the 
energy range of 10-70 keV, and GSO crystal scintillators detectors covering the 70-600 keV 
energy band. The XIS consists of four CCD detectors, three of them are front illuminated 
(FI) and one is back illuminated (BI). 

We have used the data processed with Suzaku data processing pipeline ver. 2.3.12.25. The 
observations carried out on 15 January 2009 (Obs. ID-403049010; hereafter quoted as Obs-1) 
has an useful exposure of 49.7 ks and that on 28 January 2009 (Obs. ID-403050010; hereafter 
quoted as Obs-2) has an useful exposure of 44.2 ks. In both observations, the source was 



observed at 'HXD nominal' pointing position q The XIS was operated in 'standard' data 
mode in the '| window' option which gave a time resolution of 2 s. 

For extracting the XIS and HXD/PIN light curves and spectra, we used cleaned event 
data. The cleaned events of the XIS data were screened for standard criteria Q We extracted 
the XIS light curves and spectra from these cleaned events by selecting a circular region of 
3' around the source centroid. The XIS count rates for Obs-1 and Obs-2 were 44.5 c/s and 
8.5 c/s, respectively. We also extracted the background light curves and spectra by selecting 
a region of the same size 18' away from the source centroid. Response files and effective 
area files were generated by using the FTOOLS task 'xisresp'. For HXD/PIN background, 
simulated 'tuned' non X-ray background event files (NXB) corresponding to January 2009 
were used to estimate the non X-ray background, and the cosmic X-ray background was 
simulated as suggested by the instrument teanj^ applying appropriate normalizations for 
both cases. Response files of the respective observations were obtained from the Suzaku 
Guest Observer Facility^ 

Observations of a bright source like 1A 1118- 61, especially at th e peak of the outburst 



should be checked for possible effects of pileup. ISuchy et al.l (120111 ) had applied the tool 
'pile_estimate.sl' to estimate the pileup fraction in XIS and had excluded regions of the CCDs 
with pileup fraction above 5%. We checked for the effect of pileup for Obs-1 by extracting 
spectra with the central 10 " and 15 " radii removed, and generated their corresponding 
responses. Comparing these two spectra with the original spectrum, we observe that the 
spectral shape remains almost the same in all the cases and the change in photon index 



http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/abc/ website for more information on instruments 
http://heasarc.nasa.gov/docs/suzaku/analysis/abc/node9.htmyfor standard screeing criteria 



2 See 

3 See 

4 http:/ /heasarc. nasa.gov/docs/suzaku/analysis/pinbgd. html 

5 http:/ /heasarc. nasa.gov/docs/suzaku/analysis/pin_cxb. html 
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Figure 1. The top panel shows the one day binned RXTE-ASM (1.5-12 keV) light curve of 1A 1118-61, the bottom panel 
shows the one day binned light curve from the Swift-B AT all sky monitor (15-50 keV). The arrows in the top panel indicate 
the time of the Suzaku observations of 1A 1118-61. 

is only ~0.02 for the XISs, where this effect is the most pronounced for XIS-3. Hence, 
our results are not expected to vary significantly due to this effect, since our main focus 
is on variations of the spectral parameters with pulse phase that is much larger than this 
difference. Obs-2, being an order of magnitude fainter than Obs-1 does not require this 
check. 



3 TIMING ANALYSIS 
3.1 Pulse profiles 

For the timing analysis, we applied barycentric corrections to the event data files using the 
FTOOLS task 'aebarycen'. Light curves with a time resolution of 2 s and 1 s were extracted 
from the X1S (0.2-12 keV) and the HXD/PIN (10-70 keV), respectively. We searched for 
pulsations in both observations by applying pulse folding and \ 2 maximization technique. 
The pulse period was found to be 407.49±0.73 s for Obs-1 (MJD 54846) and 407.40±0.40 
s for Obs-2 (MJD 54859) respectively. We created energy resolved pulse profiles by folding 
the light curves in different energy bands with the respective pulse periods. Light curves 
from the three XISs were added together to create the pulse profiles in the 0.2-12 keV 
energy range. The pulse profiles are found to have complex energy dependent structures on 
both the days, with more than one peak in the XIS energy band, and a single peaked but 
asymmetric structure in the PIN energy band, as shown in Figures [2] & [31 In Obs-1 (Figure 
[2]), the profiles show a double peaked structure at lower energies (less than 12 keV). The 
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Figure 2. Energy dependent pulse profiles of 1A ff 18-6f using the XIS and PIN data for Obs-1. Captions inside the figure (in 
units of keV) represent the energy range for which the profiles were created. The pulse profiles in the 0.2-f 2 keV energy range 
are made from the XIS data (data from all the three XISs are added together), where as profiles in the 12-70 keV range are 
obtained from the PIN data. 



secondary peak (phase ~ 0.2) decreases in strength with energy and disappears at ~12 keV. 
At high energies (greater than 12 keV), the profile shows a single peaked structure with a 
large pulse fraction. A narrow low energy peak (less than 2 keV) is seen coincident with the 
pulse minima at higher energies (phase ~ 0.46). For Obs-2 (Figure |3]), one main peak (phase 
~ 0.5) is seen at all energy bands, which is narrow and symmetric at low energies (less than 
3 keV). The main peak gets broadened up to 12 keV, and becomes a broad asymmetric pulse 
at higher energies. A narrow and weaker peak (at phase ~ 0.2) appears at energies greater 
than 2 keV before the main peak. At higher energies (greater than 12 keV), this feature 
gradually merges with the rising part of the main peak. Another pulse component (phase 
~ 0.75) is seen after the main peak only in the energy band below 2 keV and 6-7 keV. The 
pulse profiles show a single asymmetric peaked structure in the energy ranges greater than 
12 keV with the pulse minima at phase ~ 0.8. The pulse profiles show the presence of an 
additional peak at phase ~ 0.96 in the lowest energy bin (0.3-2 keV). The complex energy 
dependence of the pulse profile implies significant changes in the energy spectrum with the 
pulse phase, and also a complex underlying spectrum. A detailed discussion is presented in 
section 5.1. 
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Figure 3. Same as Figure 2 for Obs-2. 



3.2 Power Density Spectra 

We created the Power Density Specta (PDS) using the FTOOL task ''pouispec' from the XIS 
light curves (0.2-12 keV) which was divided into segments of 4096 s with a time resolution 
2 s. The PDS obtained from all these segments were rebinned with a geometric rebinning 
factor, and were averaged to produce the final PDS. The average PDS was also normalized 
such that its integral gave the squared rms fractional variability and expected white noise 
level was subtracted. The PDS from both observations clearly shows the peak at ~0.0025 Hz 
corresponding to the spin frequency, with its harmonics visible at higher frequencies. Obs-1 
shows a quasi-periodic oscillation (QPO) feature at 0.087 ± 0.008 Hz. 

We fitted the PDS with a power-law component corresponding to the continuum by avoiding 
the peaks corresponding to the spin frequency and its harmonics. The QPO feature was fitted 
with a Gaussian function, giving a detection significance of 13a, quality factor [y / FMHM) 
of 4.2±0.4 and the rms fractional variability of 4.58±0.17 % in the 0.2-12 keV energy range 
of XIS. Obs-2 shows a broad feature at 0.026 ± 0.003 Hz. We have fitted the feature with 
a Gaussian function which gave a detection significance of 7a, quality factor of 1.70 ± 0.2 
and rms value of 7.93 ± 0.58%. The energy dependence of the QPO has been investigated 
for Obs-1. Multiple PDS were extracted for seven different energy ranges between 3 and 10 
keV, each with a width of 1 keV in the XIS energy band. The rms of the QPO feature shows 
an increase from 4.09% to 7.74% in the XIS energy range. Figure H] shows the PDS for both 
observations in the energy band of 0.2-12 keV. We note here that it is quite unusual to 
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Figure 4. Average PDS of 1A 1118-61 for both observations showing the QPO feature in the energy band of 0.2-12 keV. The 
left panel shows the PDS for Obs-1 and the right panel shows the same for Obs-2. In both the figures, the spin period at 0.0025 
Hz and its harmonics are clearly seen with the QPO at 0.087 Hz visible for Obs-1 and the broad feature at 0.026 Hz visible for 
Obs-2. 

detect QPOs with imaging spectrometer detectors. In this source we have not only detected 
the QPOs, but are also able to investigate it's energy dependence. 



4 SPECTRAL ANALYSIS 



4.1 Pulse phase averaged spectroscopy 

The continuum emission from HMXB accretion powered pulsars can be interpreted as Comp 
tonization of soft X-rays in the plasma above the neutron star surfa ce. It can be modelled 



as a power-law with an exponential cutoff or as a broke n power-law (j White. Swank fc Holt 



1983 



Mihara 



1995 



Coburn. Heindl fc Rothschild 



2002|) 



In the presence of strong magnetic fields, the continuum photons produced inside the 
hot region by thermal bremsstrahlung and Comptonization may get resonantly scattered in 
the line forming regions, producing a cyclotron resonant scattering feature (CRSF). Often 
a strong narrow Fe Ka line can also be found originating from the cold dense matter in the 
vicinity of the neutron star. 

We performed pulse phase averaged spectral analysis of 1A 1118-61 using spectra from 
the front illuminated CCDs (XISs-0 and 3), back illuminated CCD (XIS-1), and the PIN. 
Spectral fitting was performed using XSPEC vl2.6.0. The energy ranges chosen for the fits 
were 0.8-10 keV for the XISs, and 10-70 keV for the PIN spectra. We neglected the 1.75-2.23 
keV energy range in order to take out the artificial structures in the XIS spectra at around 
the Si edge and Au edge. After appropriate background subtraction we fitted the spectra 
simultaneously with all the parameters tied, except the relative instrument normalizations 
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which were kept free. The XIS spectra were rebinned by a factor of 5 from 0.8-6 keV and 
7-10 kev. However due to the presence of a narrow iron line, the spectrum was rebinned 
only by a factor of 2 in the 6-7 keV energy range. The PIN spectra were rebinned by 
a factor of 3 in the 10-70 keV range. We had initially tried to fit the continuum with a 
pa rtial covering h igh-ene rgy cutoff model (similar modeling for this source has been done 



in 



Devasia et al. 



( l2011af )) which gave a very low value of high-energy cutoff (~0-3 keV) 
for some of the phase bins while performing pulse phase resolved spectroscopy described 
later in section 4.2. High energy cutoff in this scenario is equivalent to the cutoff power- 
law model. Hence we fitted the broad band energy spectrum of 1A 1118-61 with a partial 
covering cutoff power-law model with an interstellar absorption and a gaussian line for 
an Fe Ka feature at 6.4 keV. In addition, a broad cyclotron absorption feature was also 
found in the spectra at 49 ± 0.5 keV in Obs-1 and 52^3 ? keV in Obs-2. This feature was 
fitted with the XSPEC standard model 'cyclabs' which is a cyclotron absorption line having 
a Lorentzian profile. Presence of the bro ad cyclotron feature at ~ 55 keV had already 



been reporte d from RXTE ob servations (jDoroshenko et al. 



observations (puchy et al. 



2010) and from these Suzaku 



20111 ). The analytical form of the model which was used to fit the 



spectra is as follows: 

F(E) = exp _ffCE)NH1 ((Nl + N2 e - a{E)Nu2 ) * CYABS(E) + G) 1(E) (1) 
Where CYABS is the Lorentzian cyclotron line profile given by 

CYABS(E) = eM-Dc s , lp {Wc ; d % E °>f 2 ] (2) 

(E - EcyclY + W C ycl 

where Ec yc i is the energy of the cyclotron line & Dc yc i & Wc y d are the depth and width of 
the cyclotron line. 1(E) can be expressed as 

I(E) = E- r exp(^) (3) 

Where T is the power-law photon index and (5 is the e-folding energy of the exponential roll- 
over in keV. Nl and N2 are the normalizations of the two power-laws. Nhi is the Galactic 
hydrogen column density along our line of sight, and Njj2 is the column density of the 
material that is local to the neutron star, a is the photo-electric cross-section and G is the 
gaussian line energy. Thus the covering fraction of the more absorbed power-law is given by 

ri„, _ N2 
^Ufract — N1+N2' 

Spectral fitting showed that Obs-1 was softer than Obs-2. Since the energy spectrum had 
a strong pulse phase dependence, we did not expect the pulse averaged spectra to fit well 
to a simple continuum model. The best-fit had a reduced \ 2 °f 1-79 for 1022 degrees of 
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Table 1. Best fitting spectral parameters of 1A 1118-61. Errors quoted are for 90 per cent confidence range. 



Parameter 


15 January (Obs-1) 


28 January (Obs-2) 


N H i (10 22 atoms cm -2 ) 


0.86 ±0.01 


61 + 001 
u - D1 -0.02 


'\ ix2 (10 22 atoms cm~ 2 ) 


14 - 25 -0.59 




Cvf rac t 


0.24 ±0.01 


0.73 ±0.01 


Powlndex (r) 


0.37 ±0.02 


n Q7+0 02 
u - M '-0.04 


E-folding energy (keV) 


2i.6l8:Z 


24.8ti« 


ECycl 


49.3i°; 5 5 


51.9 + si 


Dcycl 


1.70 ± 0.05 


i-oolS-i 


Wcycl 


19.2±1.1 


10.7l|| 


Iron line energy (keV) 


6.40± 0.01 


6.40±0.01 


Iron line eqwidth (eV) 


51.0 ±1.5 


41.9 ±4.0 


Flux (XIS) (0.8-10 keV) a 


2.45 ± 0.01 


0.45 ± 0.04 


Flux (PIN) (10-70 keV) b 


1.02 ± 0.01 


0.11 ± 0.04 


reduced x 2 /d.o.f 


1.79/1022 


1.21/1023 



a Flux is in units of 10 9 ergs cm 2 s 
b Flux is in units of 10 -8 ergs cm~ 2 s' 
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Figure 5. The pulse-phase averaged spectrum of 1A 1118-61 at the peak of the outburst (Obs-1; left panel) and at the decline 
(Obs-2; right panel), using XIS from 0.8-10 keV and PIN from 10—70 keV. The upper panel shows the best-fit spectra along 
with the model components. The residuals are given in the bottom panels. 



freedom and 1.21 for 1023 degrees of freedom for Obs-1 and Obs-2, respectively. Table 1 
shows the best-fit parameters for the model. Figure [5] shows the best-fit spectra along with 
the residuals for both observations. 



4.2 Pulse phase-resolved spectroscopy 

The strong energy dependence of the pulse profiles, as shown in Figures 2 & 3 implies a 
dependence of the energy spectrum on the pulse phase. To investigate this, we performed 
pulse phase-resolved spectroscopy by dividing the data into 25 phase bins, applying phase 
filtering in the FTOOLS task XSELECT. Plotting the ratio of the spectra at some interesting 
phases, like the pulse peaks and dips to the phase average spectrum can clearly bring out 
the complexity of the spectral changes with pulse phase. In Figures [6] & [7J we plot the ratio 
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of few of the pulse phase-resolved spectra with the pulse phase averaged spectrum for both 
observations. For phase-resolved spectral analysis, For phase-resolved spectral analysis, we 
used same background spectra and response matrices for the XIS and PIN data as was used 
in the case of phase averaged spectroscopy. We fitted the spectra in the same energy ranges 
and rebinned them by the same factor as in the phase averaged case for the XIS spectra. 
For the PIN spectra however, we rebinned Obs-1 by the same factor as the phase averaged 
spectra, but rebinned Obs-2 with a higher rebinning factor of 8 above 40 keV due to poor 
statistics. The value of Nhi and the iron line width were fixed to the corresponding phase 
averaged values (as given in Table 1). Due to limited statistics in the PIN energy range, 
we could not vary all the three parameters of the 'cyclabs' model simultaneously and hence 
fixed the cyclotron width to its phase averaged value. At the main dip of the pulse profiles, 
an additional peak at the lowest energy (less than 2 keV) was apparent in the spectra of 
both observations, that indicate a soft excess and an additional blackbody component was 
required to fit the spectra for those phase ranges. For Obs-1, the phase bin corresponding to 
0.40-0.44 required an additional blackbody component with a temperature kT=0A7 keV. 
Obs-2 on the other hand required a blackbody with a temperature kT=0.34 keV in the 
phase bin of 0.96-1.00. The radius of the region from which this soft-excess is emitted 
can be calculated as 



Where is the blackbody luminosity and T&& the corresponding temperature. The radius 
of the blackbody emitting region determined from Obs-1 and Obs-2 are ~ 6 km and 4 
km respectively. The reduced \ 2 f° r the phase-resolved spectroscopy were in the range of 
1.03-1.46 and 1.32-2.5 for 1024 degrees of freedom, for Obs-1 and Obs-2 respectively. The 
variation of the spectral parameters with pulse phase are shown in Figure El The results 
obtained from the pulse phase resolved spectroscopy as seen in this figure are summarized 
below: 

Obs-1 : There was an abrupt increase in the absorption column density (-/V# 2 ) m the narrow 
dip of the profile at phase ~0.1 although the covering fraction (Cvf rac t) remained small. In 
the main dip of the profile at pulse phase of ~0.5, both the covering fraction and the Nh2 
value increased. Spectral hardening was observed in both the peaks of the pulse profile. 
Obs-2 : 

The covering fraction remained nearly constant at a high value of ~ 0.8 throughout the 
pulse phase indicating a more or less symmetric distribution of material around the neutron 
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star. The value of Nh2 increased dramatically at the main dip (phase ~0.9), although the 
covering fraction decreased to ~ 0.5. In the phase range ~0.5-0.7 in Obs-2, no high-energy 
cutoff was required to fit the spectra. The reason for this is not very clear. 
For both observations, the power-law continuum parameters like the power-law index {P index) 
normalization (P norm ), and the cutoff energy (Ecut) varied significantly with the pulse phase. 

Phase resolved spectroscopy of the cyclotron line feature: 
The high energy part (greater than 40 keV) of the data of Obs-1 had enough statistics 
to investigate the variation of cyclotron energy and depth as a function of phase. Instead 
of varying all the three cyclotron parameters, we fixed the cyclotron width to the phase 
averaged value. To perform phase-resolved spectroscopy of the cyclotron feature, 25 phase 
bins were used for the XIS data as in the previous cases. Due to limited statistics in the 
PIN data, pulse phase resolved PIN spectra with phases centered at the XIS phase bins were 
generated but at thrice their widths, resulting in 25 overlapping phase bins out of which 
only 8 were independent. Figure [9] shows the variation of the cyclotron energy and depth 
with phase. There is an indication of anti-correlation or a phase shift between the energy 
and the depth of the CRSF feature. The cyclotron energy shows a variation of about 10 keV 
and has two peaks at the same phase of the pulse peaks in the XIS energy band. As is also 
evident from Figure [HI The depth of the CRSF also shows a double peaked structure with 
large variation by a factor of 3 in a narrow phase interval of about 0.1. 



5 DISCUSSION & CONCLUSIONS 



We performed a detailed timing and broad-band pulse phase-resolved spectral analysis of 
the two Suzaku observations of the HMXB transient pulsar 1A 1118-61 during its 2009 
giant outburst in its peak (Obs-1) and declining phase (Obs-2). Owing to the broad-band 
capability and high sensitivity of Suzaku, we have been able to investigate the complex energy 
dependence of the pulse profiles in more detail and in narrow energy bins compared to the 
previous works. In addition, we have detected a soft-excess peak (in the 0.3-2 keV range) at 



the main dip of the pulse profiles for both o bservations. Broad - 
R XTE observation s have been performed by 



Doroshenko et al. 



band 



(120101 ) 



Nespoli fc Reig) (120111 ) who reported the presence of a CRSF at ~55 keV. ISuchy et al. 



jecxrai anaiysi 
Devasia et al. 


s oi xne 
(2011a) 


~55 keV. ISuchv et al. 



(120111 ) performed broad-band spectral analysis of the same Suzaku observations using partial 



covering cutoff power law and compTT models and confirmed the presence of the CRSF 
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Cnci gy (keV) Energy (keV) 

Figure 6. Phase-resolved to phase averaged spectral ratios for some particular phase bins (see values in parentheses) at the 
peak of the outburst on 15 January 2009 (Obs-1). 



feature at ~55 keV, discussin g also the possibility of the presence of a first harmonic at 
Suchy et al.l (120111 ) also included pulse phase-resolved analysis in three phase 



410 keV. 



bins in the two peaks and in the minima of the pulse profiles where they reported some 
variations in the continuum and the centroid of the cyclotron line with pulse phase. We 
carried out a more detailed phase-resolved spectral analysis of the Suzaku observations in 
narrow phase bins using a partial covering cutoff power law model. Since the aim of the paper 
is to mainly investigate the variations of the spectral parameters with the pulse phase, we 
have not included GSO data in our analysis as the data quality is not suitable to perform 
phase-resolved analysis in narrow phase bins. We, therefore, cannot investigate the possible 
presence of a harmonic at ~110 keV. Variations of the CRSF parameters over the pulse 
phase are brought out very clearly in this source, and both the energy and the depth of the 
CRSF shows some interesting patterns of variation with pulse phase. 



Primary peak (0.48-0.52) 



r +++ + 



5 10 
Energy (keV) 

Third peak (0.68-0.72) 



5 10 
Energy (keV) 
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Energy (keV) 
Main dip (0.88-0.92) 



5 10 
Energy (keV) 



Figure 7. Same as Figure 6 but now for Obs-2. 



5.1 Pulse Profiles: Energy Dependence & Soft-Excess Peak 



Energy dependence of the pulse profiles have been found 



ered X-ray pulsars such as 4U 0115+63 (ITsygankov et al 



previ ously in other accretion pow 



20071 1 and Her X-l flNagaselll989f ) 



The energy dependence of the pulse profiles of the Be HMXB pulsar 1A 1118-61 during its 
2009 outburst shows a multi peaked structure below ~ 10 keV, and a single asymmetric 
pe ak thereafter. The s e feat u res have been reporte d pre viously using RXTE o bservations 



by 



Doroshenko et al 



(120101 



Devasia et al 



fl2011af ). and iNespoli fc Reid (1201 if ). Owing to 
the better energy resolution, and the broader energy coverage, especially below 3 keV, the 
Suzaku data show a more detailed and clearer picture of the energy dependence of the pulse 



profiles. 



Suchy et al 



(120 111 ) analyzed Suzaku data using much wider energy bins than ours to 
probe the energy dependence. The presence of dips in the pulse profiles can be explained due 
to the presenc e of an additional ab sorption component at that pulse phase which obscures 



the radiation ([Galloway et al.ll200ll ). Recently similar energy dependence of the pulse pro- 



files have also been found in other HMXB pulsars for example in GRO J1008-57 (INaik et al. 
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0.5 1 

Pulse Phase 

Figure 8. Variation of the spectral parameters with phase along with the pulse profile in the 0.2-12 keV range (obtained 
from XIS data; top panels), source flux in the 0.8-10 keV range (second panels from top), pulse profile in the 10-70 keV range 
(obtained from PIN data; third panels from top) and source flux in the 10—70 keV range (fourth panels from top). The left 
panels show the variation of the spectral parameters with pulse phase for Obs-1 and the right panels show the same for Obs-2. 
The units of flux and nH2 (equivalent hydrogen column density) are in 10 -09 erg cm~ 2 s~ 1 and 10 22 atoms cm -2 , respectively. 
Power-law normalization Pnorm is in units of photons keV~ 1 cm~ 2 s~ 1 at 1 keV 



201ll ). and GX 304-1 fjPevasia et al.l l2011bl ). where the presence of the dips in the pulse 
profiles could be explained due to the presence of additional absorption component at those 
phases. 

We have detected a soft excess peak in the lowest energies (0.3-2 keV range) at the main dips 
of the pulse profiles for both observations. A soft excess emission is common among accretion 



powe red pulsars that have a low line of sight absorption (IPaul et al. 



2002 



Hickox. Narayan. &: Kallman 



20041 ) . The presence of a soft excess in acc reting X-ray pu 



cessing in the thick magnetosphere shell (IMc Cray et al 



sars can be explained as by repro- 



19821 ) or in the accretion column. 



The soft-excess could also be due to the reprocessi ng of the hard X-rays by the optical 



1 r 



thick material in the inner edge of the accretion disk (lEndo et al.ll2000l ; iRamsay et al.ll2002l ). 



y 



In this work, fitting a blackbody component to the spectra of the phases corresponding to 
the soft excess peaks, we derived the radius of the soft excess emitting regions to be ~ 6 km 
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Figure 9. Variation of the cyclotron energy and depth with the pulse phase for Obs-1. Only 8 of the 25 points are independent. 
cycLenergy denotes the energy of the centroid of the CRSF feature & cycLdepth denotes the depth of the feature. 



and ~ 4 km for Obs-1 and Obs-2 respectively. The small radii in both the cases suggest that 
the soft excess is being produced in the accretion mound itself or as reprocessed emission 
from the accretion column. Its origin from the inner accretion disk is very unlikely since the 
size of the emission region is very small compared to the inner disc radius which is of a few 
hundred kms. Also, the soft excess produced by the reproces sed disc blackbody component 



is us ually seen over a broad phase range in other pulsars (jPaul et al 



2002 



Naik fc Paul 



2004 ) and not confined to a narrow phase range like in this source. 



5.2 Quasi Periodic Oscillations 

The results r elated to the QPO fea ture found at 0.087 ± 0.009 Hz is almost similar to the 



one found by 



Devasia et al. 



( 12011af ) including its energy dependence. If the broad feature at 
0.026 ± 0.003 Hz detected in Obs-2 is interpreted as being a QPO, assuming the generation 
of the QPO at the inner radius of the accretion disc, the decrease in the QPO frequency is 
consistent with the decrease in one order of magnitude of X-ray luminosity from Obs-1 to 



Obs-2 {v qpo oc L 



3/7x 



Using this argument, the inferred value of the magnetic field is also 



in agreement with the one found by the CRSF (jPevasia et al.ll2011al ). The quality factor of 
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the QPO feature in Obs-2 is however very low compared to the values found for the other 
accretion powered pulsars. 



5.3 Pulse Phase averaged Spectroscopy 



The strong and complex energy dependence of the pulse profiles bring out the requirement 
of a multicomponent model to fit the energy spectra of 1A 1118-61. Some of the continuum 
models that have been used extensively to describe the spectra successfully ar e the Negative 



fc Positive powerlaws with a common 



Makishima et al. 



1999 



Terada et al. 



Ferm i -Dirac cutoff energy ( 



2000 



Tanaka 



xponential (NP 



2006 



Naik et al. 



X) co ntinuum model ( IMihara 



1995 



20081 ) and power-law with smooth 



19861 ). The partial covering absorption model (jEndo et al. 



Mukherjee fc Paull 120041 1 which consists of two power-law continua with a common 
photon index but different absorption components also describes the spectrum of these 
sources well. We had initially tried to fit the data with a partial covering model with a 
high-energy cutoff, but the spectra in some phases had a very low value of cutoff energy 
(~ 0-3 keV) which is equivalent to a cutoff power-law model. The phase averaged data 
was thus well fitted with a partial covering cutoff power law model. The spectrum was 
softer and steeper at the declining phase of the outburst simila r to what was fou nd in some 



other HMXB transients like V 0332+53 (IMowlavi et al.ll2006l ) and A 0535+26 (ICaballero 



20091 ). A narrow Iron Ka line centered around 6.4 keV was detected in both observations 



which implied the presence of neutral or weakly ionized Iron. The possible origin of the 
iron line is the accretion disk or the circumstellar material, or the cool dense matter in the 
magnetosphere which reprocesses the emission. 

The presence of the CRSF feature in the spe ctra had been rep orted previously by 



Doroshenko et al. 



(1201(1 ). iNespoli fc Reigl (12011( 1. & ISuchv et all (120 111 ). From the spectral 
fitting with Lorentzian profiles, we found that the CRSF line centers at 49 ± 0.5 and 52^^ 
keV for Obs-1 and Obs-2 respectively. The best fitting centriods for Obs-1 and Obs-2 are 
consistent within confidence limits. The result is also consistent with a slight increase in the 
centroid energy with decreasing luminosity which is in agreement with the variation of the 
cyclotron line energy with that of the height of the accretion column. Higher the luminosity, 
lower in energy is the CRSF produced. This is because, higher luminosity implies higher 
accretion rate and a higher height for the accretion colu mn. The CRSF feature is thus pro- 



duced further away from the surface of the neutron star (IMihara et al. 



2004 



Nakajima et al. 
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20061 ). The CRSF feature is also broader in Obs-1 which may be due to the superposition 



of CRSF features produced in scattering regions at different heights of the accretion column 
which is higher in this case than in Obs-2. 



5.4 Pulse Phase resolved spectroscopy 



In addition to pulse phase averaged spectroscopy, we have performed pulse phase-resolved 
spectroscopy to understand the complex energy dependence of the pulse profiles. The column 
density of the material local to the neutron star (iV^ 2 ) of ~ 10 24 atoms cm~ 2 can explain 
the dips seen in the pulse profiles. The changes in the value of Nh2 and the covering fraction 
with pulse phase give us an idea about the properties of the plasma in the accretion stream. 
The plasma may be a clumpy structure with different values of opacity and optical depth. 
We find that the accretion column in Obs-1 represents a more clumpy structure with the 
dip at phase ~0.1 having a very high value of Nh2 (~ 10 24 atoms cm~ 2 ) but a low 
covering value (~ 0.4). However, both values remain high in the main dip. Further, the 
accretion column in Obs-2, shows a more symmetric distribution of matter with a more or 
less constant value of covering fraction throughout the pulse phase. 

The spectral hardening at the peaks of the pulse profile in Obs-1 may represent the pas- 
sage of the magnetic axis through our line of sight, and a deep and a more direct view into the 
emission region along the magnetic axis ( jPravdo et al.lll976l ). Si milar regions of puls e hard - 



ening have be en observed previously in other sources like Her X-l (IPravdo et al 



4U 115+63 (IJohnston et al. 
198lk 



1978 



Rose et al. 



1977 



1978) 



1979J), & GX 1+4 footv. Hoffman k Lewin 



Another feature that has been found is the variation of the CRSF parameters with pulse 
phase. Since the cross section and the depth of the cyclotron scattering feature is thought 
to depend significantly on the viewing angle of the accretion column, the CRSF properties 
may change with the rotation of the neutron star, hence with pulse phase. A large variation 
of ~ 10 keV in the cyclotron energy and a factor of 3 variation in the depth could either 
imply the dipole geometry being viewed at different angles at different pulse phases, or a 
more complex underlying magnetic field structure. Fitting of this observed variation of the 
cyclotron parameters with phase with a detailed and suitable theoretical model can be used 
to infer the structure of the magnetic field around the neutron star and also put constraints 
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on its other parameters like the inclination between the magnetic and rotation axis etc. This 
is however beyond the scope of the present work. 
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